The three-dimensional Ising model: A paradigm of liquid-vapor coexistence in nuclear 

multifragmentation 
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Clusters in the three-dimensional Ising model rigorously 
obey reducibility and thermal scaling up to the critical tem- 
perature. The barriers extracted from Arrhenius plots depend 
on the cluster size as B oc A" where a is a critical exponent re- 
lating the cluster size to the cluster surface. All the Arrhenius 
plots collapse into a single Fisher-like scaling function indi- 
cating liquid-vapor-like phase coexistence and the univariant 
equilibrium between percolating clusters and finite clusters. 
The compelling similarity with nuclear multifragmentation is 
discussed. 

Nuclear multifragmentation is a process occurring at 
the limits of nuclear excitation, and, as such, portrays 
an appropriate richness and complexity. While the fun- 
damental problem of dynamics vs. statistics is still de- 
bated, it appears ever more clearly that many ther- 
mal/statistical features underlie the empirical body of 
data. In particular, two features associated with the frag- 
ment multiplicities are found to be quite pervasive in all 
multifragmentation reactions. They have been named 
"reducibility" and "thermal scaling" [[j] ||| . 

Reducibility is the property of the n- fragment emission 
probability of being expressible in terms of an elementary 
one-fragment emission probability. This property can oc- 
cur only if fragments are created independently from one 
another and it coincides with stochasticity. Both bino- 
mial, and its limiting form, Poissonian reducibilities have 
been extensively documented experimentally jj]-^. 

Thermal scaling is the linear dependence of the loga- 
rithm of the one- fragment probability with 1/T (an Ar- 
rhenius plot). It indicates that the emission probability 
for fragment type i has a Boltzmann dependence 
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where B t is a barrier corresponding to the emission pro- 
cess. 

The combination of these two empirical features power- 
fully attests to a statistical mechanism of multifragmen- 
tation in general, and to liquid-vapor coexistence specif- 
ically |. 

Many statistical models have been proposed as an ex- 
planation for multifragmentation. It is our intention here 
to identify a model which, on one hand, is as simple as 
possible, yet on the other captures the essential features 
observed in the experiments. Percolation in its many 



versions has been widely used However, while be- 

ing simple, it does not lend itself to a non-trivial ther- 
mal study P]. The three-dimensional Ising model sat- 
isfies both the criteria of simplicity in its Hamiltonian 
and lends itself to a thermal treatment with nontrivial 
results. While the Ising model has been widely studied 
in terms of its continuous phase transition, the problem 
of clusterization, of paramount interest to us here, has 
received relatively little attention [p~0| — |l9|l - We will show 
that this model contains both features of reducibility and 
thermal scaling. In showing the features of thermal scal- 
ing we will demonstrate that the slopes of the Arrhenius 
plots associated with the individual masses of the frag- 
ments, or the barriers in Eq. ([!]), portray a dependence 
on the fragment mass (A) of the form B oc A a ', where 
(j is a proper critical exponent which relates the mass to 
the cluster surface. In addition, the individual Arrhenius 
plots for each fragment mass can be absorbed into a sin- 
gle scaling function identical to that of Fisher's droplet 
model pT|-po| , which defines the liquid- vapor coexistence 
line up to the critical temperature. 

The Hamiltonian of the Ising model has two terms: 
the interaction between nearest neighbor (n.n.) spins in 
a fixed lattice and the interaction between the fixed spins 
and an external applied field H ext : 



H = - J 



E 



(2) 
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where J is the strength of the spin-spin interaction. In 
this model, the two phases of the system are a liquid 
(clusters of sites with connected up spins) and a gas 
(clusters of sites with connected down spins). In the 
absence of an external field, the distinction between up 
and down spins vanishes. Thus in the absence of an ex- 
ternal field, these two phases are in coexistence below 
the critical temperature, T c . The critical temperature 
for the three-dimensional Ising model has not been de- 
termined analytically; however, Monte Carlo techniques 
have yielded a value of T c = 4.513 J/fcf, JD|. 

In the present study, Monte Carlo techniques are used 
to determine the equilibrium cluster distribution as a 
function of temperature for a simple cubic lattice with 
periodic boundary conditions. The lattice contains 50 3 
spins, which is large but may not be large enough to 
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FIG. 1. Ratio of the variance to the mean for distributions 
of clusters with A constituents. The dashed line corresponds 
to the critical temperature, k b T c /J = 4.513. 



represent an infinite system. The Swendsen-Wang algo- 
rithm | jl6| was used to determine the equilibrium spin 
configurations of the lattice for given temperatures, and 
clusters were identified using the Coniglio-Klein @ pre- 
scription. Details about the model can be found in ref- 
erence po[ |. 

To test if the fragment distributions produced by the 
model exhibit binomial or Poissonian characteristics, the 
ratio of the variance to the mean (a 2 / (tla)) of the frag- 
ment distributions (tia) for various fragment sizes (A) 
as a function of temperature was determined. For most 
clusters, the ratio is very close to one, which is the Poisso- 
nian limit, over the entire temperature range considered 
here, which extends well beyond the critical tempera- 
ture. Only for the very lightest clusters (A < 20), does 
the ratio have a discontinuity at the critical tempera- 
ture. However, the ratio is near one even for these light 
clusters at all temperatures away from the critical tem- 
perature. (See Fig. |l|.) Thus, for all but the smallest 
clusters, the distributions exhibit Poissonian reducibil- 
ity. This signifies that the clusters are produced nearly 
independently of one another, a feature also observed in 
percolation models, Fisher's droplet model and nuclear 
fragmentation Q. 

If the fragment distributions also exhibit thermal re- 
ducibility, the distributions must be of the form given in 
Eq. (0). Thus in an Arrhenius plot (a semi-log graph of 
the number of clusters of size A (tia) vs. 1/T), the distri- 
butions should be linear below the critical temperature. 

As shown in Fig. |^, this is indeed the case over a wide 
range of temperatures (0 < T < T c ) and fragment sizes. 
While we have shown distributions for clusters up to size 
A = 100, the trend continues for larger clusters, however 
statistics decrease significantly as the size of the cluster 
increases. This linearity extends over more than four or 
five orders of magnitude. It rigorously confirms the form 
of Eq. (Q) and signifies the independent thermal forma- 



tion of fragments controlled by a single size-dependent 
barrier. This feature has been amply verified in nuclear 
multifragmentation. By fitting the fragment distribu- 
tions below the critical temperature (J/fcbT > 0.2216), 
the barriers can be extracted. The barriers for each clus- 
ter size are shown in Fig. ||. 

These barriers should find their origin in the number of 
broken bonds associated with a cluster. Therefore they 
should be well described by a power-law: 



B(A) = B A° 



(3) 



The fit of the extracted barriers according to Eq. (|J) is 
B = (12.77 ±0.04^4(°- 7230±0 - 008 ) and is remarkably good 
as shown in Fig. |. The value for a (0.7230 ± 0.0004) 
is close to 2/3, the value one would expect for spherical 
clusters of closely packed spherical objects fl5||, and to 
0.64 the value expected from scaling relations Thus 
this is a new and effective method to determine the crit- 
ical exponent a, never used before to our knowledge. 

The features of reducibility and thermal scaling dis- 
cussed above can be found united in Fisher's formula for 
the cluster abundance in a vapor as a function of cluster 
size and of temperature. The formula is 

,ms A-T /^ A Mn ,c Q A rT . . cqA' 7 . fA . 

n A {T) = q A T exp( — ^)exp(^— - )exp(--^— ) (4) 



T 



T 



where one can see the thermal scaling up to T c and the 
dependence of the "barrier" on the cluster mass through 
the critical exponent a [^l],^2). It is interesting to ex- 
plore further the applicability of this formula to the Ising 
model. 

In addition to the linear behavior of the Arrhenius 
plots below the critical temperature, the Fisher droplet 
model also predicts that the cluster size distribution at 
the critical point must follow a power law 



7i a (T c ) = q Q A T 
where go is fixed by the normalization relationship 
Ea=i^(T c ) 
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(5) 
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Away from the critical temperature, the cluster distri- 
bution should not follow a pure power law. Thus, to 
determine r without a prior knowledge of the critical 
temperature, linear fits to the cluster distributions were 
determined for all temperatures. At the critical tem- 
perature, the fit should have the lowest \ 2 /DoF and 
thus this fit should determine r. The results for the 
best fit are shown in the lower panel of Fig. |^. From 
this method, the critical temperature was found to be 
(4.515 ± .010)J/fcs with a best fit line of the form 
n A (T c ) = (30084 ± 7)A(- 2 - 39671±0 00017 '. This temper- 
ature is consistent with the value determined for infi- 
nite systems (4.513J/fcR) and the value of r (2.39671 ± 
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FIG. 2. Arrhenius plots of the cluster distributions. Sta- 
tistical error bars are shown when they exceed the size of the 
data point. The lines are fits of the form given in Eq. (|l|). 
The critical temperature is indicated by the dashed line. 



0.00017) is close to the expected value for an infinite sys- 
tem (t°° = 2.21). The fit also determines qo which can 
be compared to the prediction of Eq. (JgJ) . The fit pre- 
dicts q£ u = 30084 while gg° = 27566. The power law fit 
of the fragment abundances at the critical temperature 
is shown in Fig. S. 

In the coexistence region, which in the Ising model 
should prevail for T < T c when the magnetization is not 
constrained, the chemical potentials of the liquid and gas 
phases are equal, thus Eq. (H) can be rewritten as: 



n A (T)A T /q - exp(-c A ff e/T) 



(7) 



with e = (T c - T)/T c . Therefore, a graph of the 
scaled cluster distributions (riA(T)A T /qo) as a function of 
eA a /T should make the distributions of all cluster sizes 
collapse onto a single curve. This scaling behavior can 
clearly be seen in Fig. ||. This nearly perfect collapse 
below the critical temperature extends over six orders of 
magnitude for the broadest range of cluster sizes and it 
is perfectly linear. Therefore the three-dimensional Ising 
model and fluids belong to the same class of universal- 
ity and can be described by Fisher's droplet model. The 
Ising clusters constructed here can be properly thought 
of as "vapor" in equilibrium with the "liquid" percolat- 
ing cluster. Coexistence of the two phases is determined 
by the observation that the empirical scaling requires 
Afj, = 0. The fact that both the three-dimensional Ising 




FIG. 3. The upper panel shows the extracted barriers from 
the fits to the cluster distributions. The line is a fit of the 
form given in Eq. (^). The lower panel shows the power law 
behavior of the cluster distribution at ktT/J = 4.515. The 
line is a fit of the form given in Eq. In both panels, error 
bars do not exceed the size of the data point. 



model and the experimental nuclear multifragmentation 
data obey the same scaling predicted by Fisher's droplet 
model indicates that nuclear multifragmentation can in- 
deed be identified as the clustering (non-ideality) in a 
nuclear vapor in equilibrium with the nuclear liquid H . 

We note that the value of a determined in this work 
is larger than 0.64, the usual value of a from studies 
of the thermal properties of three-dimensional Ising sys- 
tems. There are several possible explanations for the dif- 
ferences between the value of a determined in this study 
and the values found from other studies. First, while sta- 
tistical errors are small, systematic errors may not be so 
small. Finite size effects of the 50 3 lattice could lead to 
an overabundance of small clusters. This is caused by the 
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FIG. 4. Scaling behavior of cluster distributions. 



fact that large clusters near the edge of the lattice are un- 
able to grow beyond the lattice edges easily, even though 
periodic boundary conditions are used in this model. The 
effects of lattice size as well as other systematic errors are 
discussed in greater detail in p(J . An over-abundance of 
small clusters, with a corresponding suppression of large 
clusters, would increase the slope of the cluster distribu- 
tion leading to larger than expected values of r and qo, 
just as we have observed. 

In conclusion, we have shown that the clusterization 
in the Ising model, like nuclear multifragmentation, por- 
trays reducibility and thermal scaling. In addition, the 
Arrhenius plots allow for the extraction of "barriers" 
which are found to have a dependance of B oc A a , 
where a (which is close to 2/3) is a critical exponent. 
The reducibility and thermal scaling features in the Ising 
model can be incorporated into a Fisher-like scaling with 
A/i = 0, which is obeyed rigorously over the explored 
temperature range below the critical temperature. Thus 
the observed clusters can be interpreted as the non- 
idealities of a vapor in equilibrium with a liquid. Finally, 
nuclear multifragmentation, which is seen to share all the 
scaling observed here, should be similarly interpreted as 
the clusterization of a nuclear vapor in equilibrium with 
its liquid. 
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